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Narison's r-decay-like sum rule for determining the strange quark mass is re-investigated, tak- 
ing into account isospin-breaking corrections in the extraction of the input spectral functions from 
e + e~ — > hadrons data. The corrections, estimated using experimental data on vector meson elec- 
tromagnetic decay constants and a QCD sum rule analysis of the 38 vector current correlator, are 
shown to be especially large for the isoscalar case. The reason such large corrections are natural is 
also explained. Due to the high degree of cancellation in the original sum rule, the effect of these 
corrections on the determination of m s is significant. A new central value m s = 113 — 138 MeV is 
found, in the MS scheme at 1 GeV 2 , with significant (asymmetric) errors associated with errors in 



the input experimental data. 
' I. INTRODUCTION 

> 

\ The usually quoted values for the light (u, d, s) current quark masses are obtained by a sum rule analysis of the 
■ correlator of either (for m u + rrid) the product of two divergences of the isovector axial vector current [1,2] or (for 
m s ) the product of two divergences of the strangeness-changing vector current [3,4]. One of the potential problems 
with this approach is that, in both cases, continuum contributions to the hadronic side of the sum rule are large, 
but the continuum part of the hadronic spectral function is not known experimentally. The extracted quark masses 
thus depend crucially for their reliability on that of the theoretical ansatze for the unmeasured continuum spectral 
functions, which are constructed by analogy with an extreme form of the vector meson dominance (VMD) treatment 
Oh! of the vector isovector channel wherein the continuum spectral function is modelled as a sum of Breit-Wigner terms 
whose overall normalization is adjusted to produce the desired value at continuum threshold. The threshold value is 
\ estimated, for the pseudoscalar channel, using tree-level chiral perturbation theory (ChPT) [1], and for the S = — 1 
• scalar channel, by extrapolating data using the Omnes representation, together with experimental data on the 
Ktt scattering phases [3] . Since the relevant thresholds are many resonance widths away from the poles, the assumed 
q 2 -dependence of the resonance widths is also a crucial input. Taking the m s analysis to be specific, the "standard" 
form of the g 2 -dependent s-wave width is employed. This form results from assuming the scalar couplings of the 
resonances in question to Kir are g 2 -independent, a somewhat dangerous assumption in the scalar channel. 

The above assumptions have generally been considered plausible because an analogous version is known to allow 
a successful description of e + e~ — > 7r + 7r~ cross-sections. The analogy is, however, potentially dangerous [5-7]. In 
Ref. [7], for example, the Kix portion of the scalar, S = — 1 spectral function is obtained for all s using the Omnes 
representation of the timelike scalar Kn form factor, with K e % and Kir phase data as input. The resulting spectral 
function rises much faster just above threshold, and reaches a much lower (by a factor of ~ 3) i-Q$(1430) peak height, 
than does the model version, obtained using an assumed resonant spectral shape normalized at threshold, employed 
in Refs. [3,4]. Above s ~ 2 GeV 2 , where Ktt states no longer dominate the spectral function, a full determination 
of the relevant spectral function via the method of Ref. [7] is no longer possible; below s ~ 2 GeV 2 , however, this 
result clearly demonstrates the existence of potentially large uncertainties in the "resonance-saturation/threshold- 
normalization" ansatze for the unmeasured continuum spectral functions. This suggests it is preferable to employ a 
sum rule involving experimentally-determinable spectral data. The sum rule for m s proposed by Narison [8] is of this 
type. 
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Narison's idea is to consider the difference of isovector and hypercharge vector current correlators, II 33 — II 88 , where 
3 and 8 are SU(3) F labels, and U aa is defined by 

i J d 4 x < 0|T (j;(x)J?(0)) |0 >= (<,„<,„ - q 2 g^)U aa (q 2 ). (1) 

In the operator product expansion (OPE) for this difference, each term necessarily involves the flavor-breaking param- 
eter m s — m (where m = (m u + m d )/2). Narison proposes integrating the corresponding spectral function, weighted 
as for t decay kinematics, from threshold q 2 — 4m 2 up to a mock r mass, q 2 = m 2 . As for inclusive hadronic t 
decays [9-14] , analyticity properties allow one to re- write this integral as a correspondingly weighted integral of the 
correlator difference over a circular contour of radius m 2 in the complex q 2 plane. For large enough m t , this alternate 
representation simultaneously suppresses both contributions from the region of the contour near the positive real axis 
where perturbative QCD (pQCD) becomes unreliable, and non-perturbative contributions in the OPE relative to per- 
turbative ones [9-14]. In addition, if one ignores isospin breaking, the relevant spectral functions are experimentally 
determinable in e+e~ — > hadrons. Narison's original treatment [8], which neglected isospin breaking, produced m s 
values compatible with those of the conventional analyses [3,4] (197± 29 McV in the MS scheme, at a scale of 1 GeV 2 , 
c.f. 205 ± 19 McV [4]). Since, however, there is a high degree of cancellation between the 33 and 88 contributions 
to the sum rule, it is important to consider the possibility of isospin-breaking corrections to the hadronic spectral 
functions. We investigate this question in the present paper. 



II. NARISON'S r-DECAY-LIKE SUM RULE FOR m 3 

Narison's method is based on an analogy with analyses of the inclusive r decay ratio 

/=1 _ r 7 =i (t -> v T hadrons{~f)) 

r(T^ T p e ( 7 )) ' [> 

(for concreteness, we consider r decays mediated by the charged weak isovector vector current). R^ =1 is given by 
an integral over the J = 0, 1 scalar spectral functions of the vector isovector correlator, weighted by the appropriate 
kinematic factors [9-14], which integral can be converted into one involving the scalar correlators themselves, with 
the same kinematic weights, over the counterclockwise oriented circular contour of radius m 2 in the complex s = q 2 
plane. Narison considers the difference i? 33 — i? 88 , where i?° a (a = 3, 8) results from replacing m T by a variable mass, 
mt, and the isovector current in R ! T =1 by J£. This difference can be expressed in either the "hadronic" or "contour 
integral" representations 



[Rr\ had - ^ w f 1 * (i - (i + 10 rts) 



(3) 



[RtXontour - 67U \V ud \ 2 S EW [ ^ (l - ' ) (l + D»» , (4) 

J\s\=m 2 t m t \ m t / \ m t / 



where V u d is the ud CKM matrix element, Sew — 1-0194 is the sum of the leading-log electroweak r decay corrections 

\contour 



[15], p 33 (s) and p S8 (s) are the isovector and isoscalar spectral functions and, for sufficiently large mt, [Rt a ] C ontour can 



be evaluated using the OPE. The sum rule for m s results from equating Eqs. (3) and (4) [8]. 
For the OPE side one has, with D labelling operator dimension [8] , 



[R?-R?] OPE = \V ud \ 2 S EW ]T (S£>-W). (5) 



0=2,4, • 



The D — 2 terms result from the mass-dependent perturbative contributions to the correlators which, for a flavor- 
diagonal vector current of flavor i, are [16]: 



a (mass), \ _ 3 mi(Q ) 



8 / 17981 62 A/o , 1045 , /r , > ,^ 2x2 



+T^ a{Q2)2 (32 24C(3)) £ + ° (fl3) (6) 
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with a(Q 2 ) — a s (Q 2 )/ir and rrij(Q 2 ) the running coupling and running mass of quark j, both at scale p 2 = Q 2 = —s, 
in the MS scheme. Expanding a(Q 2 ) in terms of a(m 2 ) = a and fhj(Q 2 ) in terms of a and fhj(m 2 ) = fhj, and 
performing the resulting elementary logarithmic integrals, one obtains 



€J - ffl = 12 ^ ™ S) (1 + ha + f 2 a 2 + f 3 a 3 + • • •) 

1 1 If 



(7) 



where /i = 13/3, f 2 = 30.5846, and /„, n > 2 are unknown. Here /i differs from that of Narison (/i = 11/3), but is 
in agreement with the equal mass case of Eq. (3.8) of Ref. [11]. Narison does not quote a value for f 2 . Similarly, for 
the D — 4, 6 contributions, one finds [11,8], with p representing the deviation of the light 4-quark condensate from its 
vacuum saturation value, 

- 6 if = 367r 4 a [< m ^ s > - < m uuu >] + (mi - m 4 s ) (8) 



n 



S%-Si?= U lf 3pa if [<nu> 2 -<- ss > 2 ] . (9) 

Note that the perturbative series in Eq. (7) converges much more slowly than does the analogous mass-independent 
perturbative contribution to inclusive r decay, Pinci(a) = 1 + a + 5.2023a 2 + 26.366a 3 + • • •. For m t — m T , e.g., 
using a s (m 2 ) = 0.351 ± 0.016 [22], P incl (a(m 2 )) = 1 + 0.1114 + 0.0645 + 0.0365 + • • •, while P mass (a(m 2 )) = 
1 + 0.4828 + 0.3981 + ■ ■ -. Assumptions about the convergence of P maS s(a) based by analogy on the behavior of 
Pinci(d) can, thus, not be expected to be reliable. 

For the hadronic side of the sum rule, neglecting isospin breaking, one has, using the narrow width approximation 
for the isoscalar contributions [8], 

R? = 3 '^' 2 2 ^ f ? ds ( 1 - ^ V ( 1 + 4^ h d (10) 
27ra| M J \ m 2 J \ m 2 J m 2 e+e ^ hadr °™ 



l^\Vud\ 2 S EW 



"■EM 



i 1+ ^rj — m 2 — +(--^)+-' 



(11) 



where + • • • refers to continuum and higher resonance contributions, which are small for nit less than ~ 1.6 GeV, and 
have been estimated by Narison [8]. Unfortunately, there turns out, numerically, to be a high degree of cancellation 
(to the 10 — 15% level) in Rf 3 — i?f 8 . With m t = 1.4, 1.6 GeV, and Narison's evaluation of the hadronic side, for 
example, 

[R 33 - i?f] TO(=1 4 GcV = (1.853 ± 0.072) - (1.581 ± 0.066) = 0.272 ± 0.098 (12) 
[R 33 ~ Rt 8 ] mt=16GcV = (1.793 ± 0.070) - (1.626 ±0.069) = 0.167 ± 0.098 . (13) 

Isospin breaking at the few % level in the individual terms is, therefore, not necessarily negligible in the difference, 
especially if, as one would expect, e.g., for p-co mixing effects, the signs of the effect were to be opposite in the two 
cases. 

We now describe the input to the present version of analysis of the Narison sum rule. Modifications to the 
OPE side of the sum rule are minor. First, we use Narison's updated value of the D — 6 4-quark condensate, 
pa s < uu > 2 — (5.8 ± 0.9) x 10~ 4 GeV 6 [18], and, in evaluating the difference of light and strange quark terms, allow 
< ss > / < uu > to vary between 0.7 and 1, as in Ref. [7]. Second, for the (very small) D = 4 terms, we take the 
light quark condensate to be given by the GMOR relation, and the strange quark condensate by the kaon version 
thereof, multiplied by a factor, ck, (varied between 0.5 and 1 to account for possible SU(3)f breaking). The ratio 
m s /m = 24. 4± 1.5, determined from ChPT to 1-loop [19], is also used. Third, since the 4-loop (3(a) [20] and 7(a) [21] 
functions are now available, we employ throughout the 4-loop expressions for a(Q 2 ) and fh(Q 2 ), fixing the 3- flavor scale 
parameter, A 3 from the recent ALEPH r decay data analysis [22]. The slow convergence noted above for P mass (a) 
at m t = m T is even more accentuated at lower scales (e.g., for m t — 1.2, 1.4 GeV, P m ass = 1 + 0.672 + 0.736 + • • • 
and 1 + 0.580 + 0.548 + • • •, respectively). An estimate of the 0(a 3 ) term in P mass (a) thus appears crucial. We 
estimate using the procedure of Ref. [23] (CKS). For the three cases where the 0(a 3 ) coefficient of quadratically- 
mass-dependent observables are known, the resulting estimates are accurate to ±25 [23]. The two possible versions, 
labelled FAC and PMS, yield [f 3 ] FAC = 288.0 and [/ 3 ] PM5 = 290.1. In contrast, [1,1] and [0,2] Pade estimates yield 
[/ 3 ] [M1 = 215.9 and [/ 3 ][°' 2 l = 183.7. To be conservative, we take / 3 = 200 ± 200. 
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On the hadronic side the major change is that we now make isospin-breaking corrections to the extracted isovector 
and isoscalar spectral functions. These turn out to be quite large. Such corrections are unavoidable in the isoscalar 



where e + e 



hadrons is the only source of experimental data. In the isovector case, we will combine e + e 



hadrons and r decay data in order to reduce the errors on Rf d . For more details on the isospin breaking corrections 
to the e + e~ — ► hadrons data, see Ref. [24]. The necessity of such corrections is obvious. Just as isospin breaking is 
observed through p-ui interference in e + e~ — > hadrons, so the vector meson electromagnetic (EM) decay constants, 
Fy M , will have both isospin-conserving and isospin-violating pieces. Defining F v by < 0|J°|V, A >= F v mve\ for 
a = 3, 8, we see that Fy M = F v + -y^F v . Ff, F 3 and F? vanish in the isospin limit, but will be non-zero in the real 



(3)n2 



world. The p contribution to the EM spectral function then contains an isospin-conserving piece proportional to (F f 
and an isospin-violating piece proportional to Fp 3 ^ Fp 8 \ The latter (associated with the 38 portion of p EM ) must be 
excluded in determining the p contribution to p 33 . Similarly, contributions proportional to F V F V (V = <f>, u) should be 
removed from the physical EM widths to obtain the u> and <j> contributions to the 88 spectral function. The evaluation 
of the (unmeasured) isospin-breaking decay constants is accomplished by performing a QCD sum rule analysis of the 
mixed-isospin correlator < 0|T (J 3 J*) |0 > [24,26], for which the resonance contributions are directly proportional to 
the product F V F V . Combining the results for these products with the experimental values for Fy M , one extracts F v 
and Fy separately. The sum rule provides good constraints on the product for the p and u>, weak constraints for the <f>, 
and is insensitive to higher resonance contributions. (Details of the analysis, in relation to CVC tests and the extraction 

of the sixth order ChPT low-energy constant, Q [27,28], may be found in Ref. [24].) One finds i^ 8) = 2.8 ± 1.1 MeV, 



r(3) 



4.2 ± 1.5 MeV and f£ ] = 0.21 ± 0.21 MeV ( c.f., F 



EM 
P 



154 ± 3.6 MeV, F EM = 45.9 ± 0.8 MeV and 
F EM = —79.1 ± 2.3 MeV) [24]. These results satisfy several physical naturalness criteria [24]. The p, uj and 
contributions to the 33 and 88 vector current spectral functions are given by [p aa (9 2 )]y = [F v ] 2 <5y(s) where <V(s), 
V = u>, <f>, in the narrow width approximation, is the usual S function, while S p (s) is the corresponding p Breit-Wigner. 
The standard extractions, in contrast, are obtained by replacing F 3 with F EM and F v with ^/3Fy M for V = u, <j>. 
The corrections necessary to produce the true resonance contributions to the 33 and 88 spectral functions are thus 




0.979 ±0.0086 



1.189 ±0.065 



= 1.0054 ± 0.0054 



(14) 



The size of the overestimate in the case of the 33 spectral function is still noticeably smaller than the ~ 5% errors 
on the e+e~ — > hadrons cross-sections in the resonance region. Note that the scale (~ 1%) of the isospin-breaking 
contribution to F EM corresponds, as one might expect (since the assumption of mixing dominance corresponds to 
a leading chiral order approximation [24]), to what is obtained by assuming dominance by p-u mixing, and then 



evaluating this mixing following the updated analyses of e + e 



discussed in Refs. [251. Note also that, in 



this approximation, and neglecting both the p width and the difference of the p and ui masses, the p and u correction 
terms would cancel in spectral integrals. The approximate cancellation of the corrections associated with the results 
of Eq. (14), when one forms a sum rule involving the sum of p and lo contributions, is thus simply a manifestation 
of the dominant role of p-oj mixing. Without understanding this point, it is easy to be misled into thinking that the 
scale of the individual resonance isospin-breaking decay constants is set by that of the isospin-breaking terms in the 
corresponding OPE representation. The latter, however, is associated with the sum, which involves a rather close 
cancellation, rather than the scale of the individual terms. While the 1% correction in the case of the p is, as just 
explained, quite natural, the ui correction to the 88 spectral function (~ 19%) might, in contrast, seem unnaturally 
large to some readers. It is, however, rather easy to see that such a large correction is actually expected in this 
case. As is well-known, in the limit that the vector meson nonet is ideally mixed, but the octet vector current matrix 
elements are otherwise given by SU(3)f, the EM decay constant of the 1=1 component of the p, F£, is 3 times that 
of the 1 = component of the u, F^. Writing p = pj + ew; and u) = wj — epi, (I denoting the isospin pure states), 
with e ~ 0(6m), the physical EM decay constants become 



pEM — pi i £ pl 



7 EM 



F 1 



K (1 - 3e) . 



(15) 
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The fractional correction for the u> should thus be ~ 9 times that for the p, and of opposite sign. Both features are 
present in the results of the sum rule analysis. 

In reanalyzing the sum rule, we would, ideally prefer to work with m t near m T where the perturbative series 
Pmass{a) is better behaved. Unfortunately, at such scales, the oj' contribution to R 88 becomes important. Since the 
sum rule employed to estimate the isospin-breaking decay constants is insensitive to the p'-uj' region, we arc unable 
to correct the oj' contribution for isospin breaking, and must thus work at scales below m t ~ 1.6 GeV, where this 
contribution is small (< 2% [8]). The slow convergence of P mass {a), similarly forces us to values of m t above <~ 1.4 
GeV. 

The corrected values of R 33 , R 88 , their difference, and the resulting values of m s , are given, as a function of mt, 
in Table 1. R 33 is obtained by combining information from e + e~ — > hadrons with that from t decay. In the latter 
case, we have used the recent ALEPH tabulation of the isovector vector spectral function [29], re- fitting the portion 
of the unfolded distribution from threshold up to s = (1.6 GeV) 2 relevant to our analysis. The results turn out to be 
almost identical to those obtained by direct numerical integration of the unfolded ALEPH distribution. Although the 
errors associated with the extraction of the isovector spectral function from r decay data are a factor of <~ 2 smaller 
than for e + e~ — > hadrons, one must bear in mind that there remain at present unknown 0{olem) corrections to the 
relation between the spectral function extracted in r decay and that appearing in e+e~ — > hadrons [30]. Marciano 
[30] has assigned an additional uncertainty of ~ 3% to account for these corrections. We have added this error in 
quadrature with that quoted by ALEPH. For the estimate based on e+e~ — ► hadrons data, we employ Narison's 
evaluation of the uncorrected hadronic integral for R 33 , and PDG96 [31] values for the ui and <j> EM widths to obtain 
the uncorrected version of R 88 . Quoted values for m s in Table 1 are in the MS scheme, at scale \i = 1 GeV, using 
the 4-loop running, and correspond to central values for all input. The errors quoted for R 88 reflect both those from 
the PDG96 partial widths and those from the uncertainties in the theoretical analysis of the mixed-isospin correlator 
sum rule. As expected, the isospin-breaking corrections have a very significant impact on R 33 — R 88 - The decrease in 
the value of the difference also magnifies the effect of the experimental errors. The errors quoted for m s correspond 
to those on the hadronic integrals only (the first to that on R 33 , the second to that on Rf 8 ). Errors associated with 
uncertainties in the remaining inputs are as follows: (1) for ck- ~ ±0.4 MeV; (2) for the light-flavor four-quark 
condensate: ~ ±4 MeV; (3) for < ss > / < uu >: < ±17 MeV; (4) for / 3 : < ±14 MeV; (5) for A 3 : < ±3 MeV. 

We see from Table 1 that the corrections to R 88 significantly lower the hadronic side of the sum rule, and hence 
m s . Taking an average of the three determinations, we obtain 

m s (l GeV) = 113 ±|| +H ± 23 MeV (16) 

the first two sets of (asymmetric) errors being associated with the experimental input and the last set with uncertainties 
in the input on the OPE side. The result, while significantly below that of conventional sum rule analyses, is compatible 
with that of Colangelo et al. [7], and within errors, also the low values obtained in some recent lattice analyses [32,33]. 
To understand that the errors associated with the experimental input are still very significant, note that, if we use 
only the r decay data to determine the isovector contributions, the central values of m s are changed to 148.5, 131.5 
and 134.5 MeV, for m t — 1.4, 1.5 and 1.6 GeV, respectively. We should also point out that the existing extraction of 
m s based on flavor-breaking in hadronic t decay [34] is incorrect as a result of an error in the input coefficients of the 
perturbative series given in Ref. [16] for the mass-dependent D=2 terms entering the strangeness-changing current 
contribution to these decays [35] . 

III. SUMMARY 

We have shown that isospin-breaking corrections significantly alter the value of p 33 — p 88 extracted from existing 
e + e _ — > hadrons data, and that the resulting changes to the Narison sum rule for m s produce a central value fh s ~ 110 
MeV (at a scale 1 GeV) (~ 140 MeV if one uses only r decay data for the isovector input to the sum rule) . The errors 
on this result, which are not insignificant, are currently dominated by the errors on the experimental input. Improved 
experimental data would, of course, significantly reduce the errors on the hadronic side of the sum rule. An improved 
treatment of both the hadronic and OPE sides of the mixed-isospin vector current sum rule is also highly desirable, 
both as a check of the stability of the results for, and a means of potentially reducing the errors on the determination 
of, the isospin-breaking vector meson decay constants. 
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TABLE I. Hadronic input and extracted values of m 3 



m t (GeV) 


R? 




R 33 R S8 
tit ~ tit 


m s (MeV) 


1.4 


1.890± 0.098 


1.706± 0.057 


0.184± 0.113 


i oq i +3 / +^3 
i^o.±_ 54 _ 28 


1.5 


1.842± 0.063 


1.737± 0.056 


0.105± 0.084 


1 04 4 + 36 + 33 
_ 5 g_49 


1.6 


1.800± 0.048 


1.713± 0.056 


0.087± 0.074 


110 4 + 33 + 39 
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